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Abstract 



We study Higgs boson decays in the minimal supersymmetric standard model 
where the tree-level CP invariance of the Higgs potential is explicitly broken 
by loop effects of soft CP-violating Yukawa interactions related to scalar 
quarks of the third generation. The scalar-pseudoscalar mixing among two 
neutral CP-even Higgs bosons and one CP-odd Higgs boson due to explicit 
CP violation modifies their tree-level couplings to fermions, to the and 
Z bosons and to Higgs bosons themselves significantly. We analyze the phe- 
nomenological impact of explicit CP violation on the branching ratios of the 
neutral Higgs boson decays in detail and discuss how to directly confirm the 
existence of explicit CP violation through t + t~ and tt spin correlations in 
the decays of the neutral Higgs bosons into a tau-lepton pair and a top-quark 
pair. 



PACS number(s): 14.80.Cp, 13.88.+e 
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I. INTRODUCTION 



Revealing the physical mechanism responsible for the breaking of the electroweak sym- 
metry is one of the most crucial issues in particle physics. If the fundamental particles such 
as leptons, quarks and gauge bosons remain weakly interacting up to very high energies, the 
sector in which the electroweak symmetry is broken must contain one or more fundamental 
scalar Higgs bosons with masses of the order of the symmetry-breaking scale, However, if 
the Standard Model (SM) is embedded in a Grand Unified Theory (GUT) at high energies, 
the natural scale of the electroweak symmetry breaking would be expected close to the unifi- 
cation scale. Supersymmetry |T| provides an elegant solution to this hierarchy problem and 
furthermore leads to an excellent agreement between the value of the electroweak mixing 
angle sin 2 6y/ predicted by the unification of the gauge couplings and the experimentally 
measured value. 

Of course, supersymmetry must be (softly) broken to be phenomenologically viable. In 
general, this breakdown introduces a large number of unknown parameters, many of which 
can be complex ||. CP-violating phases associated with sfermions of the first and, to a 
lesser extent, second generations are severely constrained by bounds on the electric dipole 
moments of the electron, neutron and muon. However, there have been several suggestions 
to evade these constraints without suppressing the CP-violating phases. One option 
is to make the first two generations of scalar fermions rather heavy so that one-loop EDM 
constraints are automatically evaded. As a matter of fact one can consider so-called effective 
SUSY models [[| where de-couplings of the first and second generation sfermions are invoked 
to solve the SUSY FCNC and CP problems without spoiling the naturalness condition. 
Another possibility is to arrange for partial cancellations among various contributions to 
the electron and neutron EDM's ||. 

Following the suggestions that the phases do not have to be suppressed, many impor- 
tant works on the effects due to the CP phases in the minimal supersymmetric standard 
model (MSSM) have been already reported; the effects are very significant in extracting the 
parameters in the SUSY Lagrangian from experimental data |J, estimating dark matter 
densities and scattering cross sections and Higgs boson mass limits |pH|, CP violation in 



the B and K systems [fL0 |, and so on. In particular, it has been found || that the Higgs- 
sector CP violation induced via loop corrections of soft CP-violating Yukawa interactions 
may drastically modify the couplings of the light neutral Higgs boson to the gauge bosons. 
As a result, the current experimental lower bound on the lightest Higgs boson mass may be 
dramatically relaxed up to a 60-GeV level in the presence of large CP violation in the Higgs 
sector of the MSSM. 

Since the experimental observation of scalar Higgs particles and the detailed confirmation 
of their fundamental properties are crucial for our present understanding of the mechanism 
of the electroweak symmetry breaking, a very precise prediction of the production cross 
sections and of the branching ratios for the main decay channels is mandatory. So, several 
recent works have made the detailed predictions of the Higgs mass spectrum and their 
production cross sections in the MSSM with explicit CP violation. Along with the mass 
spectrum and cross sections, it is also very important to estimate the branching ratios of the 
main Higgs boson decays precisely. This indispensable requirement leads to our systematic 
study of the Higgs boson decays in the MSSM [ll|] with explicit CP violation, where both 
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Higgs-boson masses and their couplings to fermions, gauge bosons are significantly affected 
by loop corrections of soft CP-violating Yukawa interactions to the Higgs sector. We provide 
a detailed quantitative estimate for the branching fractions of all possible two-body decay 
modes of three neutral Higgs bosons and investigate the possibility of measuring CP violation 
through the fermion spin correlations in the neutral Higgs decays into a tau-lepton pair 
and/or a top-quark pair. 

The organization of the paper is as follows. In Section |TJ we give a brief review for 
the explicit CP violation in the Higgs sector due to the significant contributions of radiative 
corrections from trilinear Yukawa couplings of the Higgs fields to scalar top and scalar bottom 
quarks. This review is mainly based on the work by Pilaftsis and Wagner ||. The parts 
involving other supersymmetric particles such as charginos and neutralinos can contribute 
to the effective Higgs potential. But, the chargino and neutralino contributions are degraded 
in the present analysis with their masses taken to be of the order of a large SUSY breaking 
scale MgusY- Then, we introduce the general mass matrix dictating the mixing among three 
neutral Higgs states and the resulting Lagrangian describing the couplings of the Higgs 
bosons to fermions, gauge bosons, and Higgs bosons themselves. In Section [TTI| we discuss 
the effects of the CP phases on (almost) all the two-body Higgs boson decays in detail, 
which will give an important phenomenological impact on Higgs-boson searches. Section 
P"V| is devoted to investigating the possibility of detecting CP-violating effects in the neutral 
Higgs boson decays into tau-lepton pair or top-quark pair in which the subsequent decays 
of tau leptons and top quarks enable one to extract the information on the spin of the 
parent fermions; tau leptons and top quarks. Finally, Section [V] summarizes our findings 
and concludes. 



II. CP VIOLATION IN THE MSSM HIGGS SECTOR 

The MSSM introduces several new parameters in the theory that are absent from the 
SM and could, in principle, possess many CP-violating phases. Specifically, the new CP 
phases may come from the following parameters: (i) the higgsino mass parameter /z, which 
involves the bilinear mixing of the two Higgs chiral superfields in the superpotential; (ii) the 
soft SUSY-breaking gaugino masses M a (a = 1, 2, 3), where the index a stands for the gauge 
groups U(l)y, SU(2)l and SU(3)c, respectively; (iii) the soft bilinear Higgs mixing masses 
mf 2J which is sometimes denoted as Bfi in the literature; (iv) the soft trilinear Yukawa 
couplings Af of the Higgs particles to scalar fermions; and (v) the flavor mixing elements 
of the sfermions mass matrices. If the universality condition is imposed on all gaugino 
masses at the unification scale Mx, the gaugino masses M a have a common phase, and if 
the diagonal boundary conditions are added to the universality condition for the sfermion 
mass matrices at the GUT scale, the flavor mixing elements of the sfermions mass matrices 
vanish and the different trilinear couplings Af are all equal, i.e. Af = A. 

The conformal-invariant part of the MSSM Lagrangian has two global U(l) symmetries; 
the U(1)q Peccei-Quinn symmetry and the U(l)# symmetry acting on the Grassmann- 
valued coordinates. As a consequence, not all CP-violating phases of the four complex 
parameters {/i, m^ 2 , M a , A} turn out to be physical, i.e. two phases may be removed by 
redefining the fields accordingly Employing the two global symmetries, one of the 
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Higgs doublets and the gaugino fields can be rephased such that M a and m\ 2 become real. 
In this g(/i) and arg(A) are the only physical CP-violating phases in the low-energy 

MSSM supplemented by universal boundary conditions at the GUT scale. Denoting the 
scalar components of the Higgs doublets Hi and H 2 by Hi = — zt 2 $* (r 2 is the usual Pauli 
matrix) and H 2 = $ 2 , the most general CP-violating Higgs potential of the MSSM can be 
conveniently described by the effective Lagrangian 

C v = + ^($^2) + m 2 i 2 ($\$ 2 ) + mlK^i) 

+Ai($I$i) 2 + A 2 ($^ 2 ) 2 + A 3 ($ t 1 $i)($ 2 $ 2 ) + A 4 ($ t 1 $ 2 )($ 2 $i) 
+A 5 ($ t 1 $ 2 ) 2 + A;($ 2 $!) 2 + A 6 ($I$i)($t$ 2 ) + Ag($l$i)($^$i) 
+A 7 ($ 2 $ 2 )($ t 1 $ 2 ) + A;($ 2 $ 2 )($ 2 $i) . (1) 

In the Born approximation, the quartic couplings Aj.,2,3,4 are solely determined by the gauge 
couplings and As^j are zero. However, beyond the Born approximation, the quartic cou- 
plings A5 ( 6 7 receive significant radiative corrections from trilinear Yukawa couplings of the 
Higgs fields to scalar-top and scalar-bottom quarks. These parameters are in general 
complex and so lead to CP violation in the Higgs sector through radiative corrections. 
The explicit form of the couplings with the radiative corrections included can be found in 
Refs. 

It is necessary to determine the ground state of the Higgs potential to obtain physical 
Higgs states and their self-interactions. To this end we introduce the linear decompositions 
of the Higgs fields 

$1 = ( ft* A + iai) ) ' $2 = e * ( ^(t* A + ia 2 ) ) • ( 2 ) 

with Vi and v 2 the moduli of the vacuum expectation values (VEVs) of the Higgs doublets 
and £ their relative induced CP-violating phase. These VEVs and the relative phase can 
be determined by the minimization conditions on £y, which can be efficiently performed 
by the so-called tadpole renormalization techniques 0. It is always guaranteed that one 
combination of the CP-odd Higgs fields ai and a 2 (G° = cos (3ai — sin (3a 2 ) defines a flat 
direction in the Higgs potential and so it is absorbed as the longitudinal component of the 



Z boson. [Here, sin/? = v 2 j \jv\ + v 2 and cos/? = v\/ yt> 2 + v 2 ] As a result, there exist one 
charged Higgs state and three neutral Higgs states that are mixed in the presence of CP 
violation in the Higgs sector. Denoting the remaining CP-odd state a = sin/fox + cos j3a 2 , 
the 3x3 neutral Higgs-boson mass matrix describing the mixing between CP-even and 
CP-odd fields can be decomposed into four parts in the weak basis (a, <pi, <p 2 ) : 



2 _ . M 2 P M 2 PS 1 (3) 
M ° ~ I M 2 S p M 2 1 ' 1 j 



where Ml and M 2 S describe the CP-preserving transitions a — * a and (<fii,<fi 2 ) — > (^1,^2)5 
respectively, and M 2 PS = (M 2 SP ) T contains the CP-violating mixings a <-> (<pi,<p 2 ). The 
analytic form of the sub-matrices is given by 



4 



M 2 P = 



m l = 



M 2 



SP 



2 (l(X 5 e 2 ^)sp + l(X 6 e^)cp 
~ V {l(X 5 e^)cp + I(X 7 e^)sp 



2n(\ 5 e 2 *) spCp + l -K{\^)cl + ^(A 7 e 4 «)s 



Ml = 




)s 2 + 2K(\ 6 e*)s c f} 



A34W + ft(A 6 e*)cj| + ft(A 7 e* 



A34 W + ft(A 6 e*)(-| + ^(A 7 e^)4 
2X^1 + 2^(A 5 e 2i «)4 + 2H(X 7 e^)s fi c fi 



where cp = cos /3 and sp = sin /? and the mass squared m 2 is given by 



TZ(ml 2 e 



) + v 2 



2K(X 5 e 2i Z)spCp + -n(X 6 e^)cl + -TZ{X 7 S)s) 



Correspondingly, the charged Higgs-boson mass is given by 



(4) 



(5) 



m 2 H± 



TZ{X 5 e 



2 * V . 



(6) 



Taking this very last relation between m#± and m a into account, we can express the neutral 
Higgs-boson masses as functions of m H ± , fx, A t , A b , a common SUSY scale M SU sy, tan (3 and 
the physical phase £. However, with the chargino and neutralino contributions neglected, the 
radiatively induced phase £ can be absorbed into the definition of the fi parameter. Clearly, 
the CP-even and CP-odd states mix unless all of the imaginary parts of the parameters 
A5,6,7 vanish. Since the Higgs — boson mass matrix j\A\ describing the scalar-pseudoscalar 
mixing is symmetric, we can diagonalize it by means of an orthogonal rotation 0\ T = 
diag(m# , m# , m# ) with the ordering of masses m Hl < mH 2 < rn Hs . The neutral Higgs- 
boson mixing affects the couplings of the Higgs fields to fermions, gauge bosons, and Higgs 
fields themselves as shown in the following. On the other hand, there could exist CP- 



violating vertex correctios, but they have been shown to be rather small |14[ so that those 
effects are not included in the present work. 

Firstly, the interactions of the neutral Higgs fields with SM fermions are described by 
the Lagrangian 
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(7) 
(8) 



Obviously, the Higgs-fermion-fermion couplings are significant for the third-generation 
fermions, t, b and r because of their relatively large Yukawa couplings. We readily see 
that the effect of CP-violating Higgs mixing is to induce a simultaneous coupling of Hi 
(i = 1,2,3) to CP-even and CP-odd fermion bilinears // and fijsf |[L5|| . This can lead 
to a sizable phenomenon of CP violation in the Higgs decays into polarized top-quark or 
tau-lepton pairs |16|], which will constitute the topic of Section [TV[ 

Secondly, the couplings of the Higgs fields to W and Z bosons may be read off by the 
Lagrangian 
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£hvv = gmw(c(i0 2 ,A-i + spOz A -i) H { 



1 



2 Z v Z ^ 



C-hh±ht = ~{c p 3 ^i - s p 2A -i + zOi,4_i) T^ +At ( flii d„ H~ ) + h.c. 
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HHZ 



Ac w 



Oi^-ii^CpOs^-j — SgOi^-j) — 0\^-j{cp0 3) A-i — S/301,4- 



(9) 
(10) 
(11) 



Note that the Z boson can only couple to two different Higgs particles. The reason is that 
Bose symmetry forbids any antisymmetric derivative coupling of a vector particle to two 
identical real scalar fields. On the other hand, the orthogonality of the mixing matrix O 
and the assumption detO = 1, which we can take without loss of any generality, leads to 
the following important relation between the couplings of the neutral Higgs bosons to the 
gauge bosons: 



9k ^ijkQijy 



(12) 



with the definitions; g { = c p 2i + sp0 3i and = Ou(cp0 3j - sp0 2j ) - Xj {cp0 3i - sp0 2i ). 
One intermediate result from the relation and the unitarity constraint is that the knowledge 
of two gi is sufficient to determine the whole set of couplings of the neutral Higgs to the gauge 
bosons fl7| . As will be seen in Sec. Ill C, the above relation leads to a close correlations 
between B{H 2 — > H X Z) and B{H 3 — ► VV) and it also leads to a similar close relation 
between B(H 3 -> VV) and B(H 2 -> H X Z). 

Thirdly, the trilinear Higgs self-couplings also are affected by the radiative corrections. 
Their relevant interactions can be read off by the Lagrangian 



C 



HHH 



AijkHiHjHk + BiHiH + H 



where the totally symmetric Auk and the coefficients Bi are given by 



a/37=l,2,3 



Bi — O a ^ib a , 

a=l,2,3 



(13) 



(14) 



where the totally-symmetric coefficients a a ^ and the coefficients bi can be directly obtained 
from taking every possible combination of third-order derivatives of the general Higgs po- 
tential ([!]) with respect to the fields {a, 0i, <p 2l H ± }. The coefficients, of which the form is 
presented in the Appendix, are functions of the quartic couplings Aj (i = 1 to 7). As can 
be checked in the Lagrangian, the radiative corrections generate various new interaction 
vertices among the neutral Higgs fields and the charged Higgs fields, which are absent in the 
CP-invariant MSSM. 

Finally, we emphasize that the size of CP violation due to radiative corrections to the 
Higgs potential is characterized by a dimensionless parameter rjcp 



m, 



( 






\ 


v327r 2 M s 2 USY y 



sin $ 



(15) 



where $ = arg(Affj,)+^, i.e. the sum of three CP-violating phases. So, for \fi\ and \Af \ values 
larger than the SUSY-breaking scale M SU sy, the CP-violating effects can be significant. 
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III. HIGGS BOSON DECAYS 



A. Higgs boson masses 

Higgs-boson masses are not determined a priori within the theory and their decay pat- 
terns depend strongly on the masses as well as several SUSY parameters. So, in order 
to determine the possible decay modes and branching ratios, it is necessary to investigate 
the change of the mass spectrum with respect to the SUSY parameters. As explained in 
Section 0, the mass spectrum is determined by nine real quantities tan/5, m#±, \A t \, 
| Aft |, $ M , $A t , &A b , and Msusy- For simplicity, we assume in our numerical analysis for the 
following dimensionful parameters 

| At | = \A b \ = ITeV, |/i| = 2TeV, M SUS y = 0.5 TeV , (16) 

while we treat m#± as a dimensionful free parameter. Noticing that the effects of the CP 
phases appear with the unique combination $ = arg(Afn) +£ and using the freedom of field 
redefinitions guaranteed only if chargino and neutralino contributions are neglected, we take 
for the phases 

% + Z = 0, ®A t = $A b = $- (17) 

Then, we are left with two dimensionless free parameters tan/3 and $. So, eventually, 
all the mass spectrums and couplings are determined completely by three parameters 
{tan j3, $, m H ±}. 

Figure 1 shows the masses of three neutral Higgs bosons and one charged Higgs boson 
with respect to the charged Higgs mass m#± for different values of the CP phase $ with the 
fixed values for other SUSY parameters in Eq. flTE|). The value of tan (3 |T5| is taken to be 3 



and an artificial experimental condition > 70 GeV is imposed. The lower solid line is 
for the lightest Higgs boson Hi, the upper solid line for the heaviest Higgs boson H%, and 
the dashed line for the intermediate Higgs boson H2. Comparing the mass spectrums of the 
whole frames, we find the following interesting features for the Higgs boson masses: 

• The mass of the lightest Higgs boson Hi is always less than 130 GeV irrespective of 
the value of $. 

• The intermediate Higgs-boson mass is very sensitive to the CP-violating phase $ and 
becomes very close to the lightest Higgs boson mass around $ ~ 70°. 

• The heaviest Higgs-boson is almost degenerate with the charged Higgs boson except 
the region with a small charged Higgs boson mass. 

• Except the lightest Higgs boson mass, the other Higgs boson masses increase with the 
charged Higgs mass, which is treated as a free parameter. 

One can also notice that the the minimum values of the allowed heavier Higgs masses under 
the constraint niu 1 > 70 GeV depend strongly on the CP phase. These features clearly 
suggest that while the decay pattern of the lightest neutral Higgs boson might be rather 
insensitive to the CP phase and the charged Higgs mass, the decay patterns of the other 
neutral Higgs bosons and the charged Higgs boson itself are very sensitive to them. 
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B. Higgs decay branching ratios 



In most cases, the most important decay channels of the Higgs bosons are two-body 
decays to the heaviest particles because the Higgs couplings to the SM particles are propor- 
tional to their masses. Sometimes, because of the coupling enhancement below-threshold 
three-body decays of the Higgs bosons may also be important []B5|. Nevertheless, in this 
work we concentrate on only two-body decay channels. On the other hand, if supersym- 
metric particles are light enough, Higgs bosons can decay into supersymmetric particles. 
However, since we have assumed a common SUSY scale of the order of 1 TeV, decays to 
sfermions, neutalinos, and charginos are not expected to play an important role in the Higgs 
mass range of a few hundred GeV which we are analyzing. 

Let us now present the explicit form of each Higgs boson decay width. First of all, the 
decay width for the Higgs boson decays into a fermion-pair is given by 



T{Hi - ff) 



N C G 



F 2 

m f m Hi 



7T 



{(3 f v)f + (IlWj 



(18) 



where N c = 3 and 1 for quarks and leptons, and (3f is the fermion velocity J 1 — Am'j/m'jj. in 
the rest frame of the Higgs boson H^. Among these fermionic decay channels, the bottom- 
quark and tau-lepton modes give most dominant contributions to the decay width of the 
lightest Higgs boson due to their relatively large Yukawa couplings. Numerically, the branch- 
ing ratio of the bottom-quark channel is almost 90% and the tau-lepton mode is 10% in 
the almost whole range of the allowed lightest Higgs boson mass. Moreover, if the masses of 
the heavier Higgs bosons are smaller than the threshold of two W's, they will be the main 
decay channels for the heavy states as well. As stated in the expression, the fermionic decay 
width is strongly dependent on the fermion mass so that the exact determination of the 
mass is required for the precise estimate of the branching ratio. Keeping this point in mind, 
we simply take m;, = 4.25 GeV and m c = 1.25 GeV for our numerical analysis. 

Secondly, the decay widths of the Higgs bosons into a gauge boson pair can be obtained 
by the Lagrangian for the Higgs couplings to the gauge bosons. The analytic form of those 
decay widths can be summarized as follows: 



T(Hi -> VV) 



G F m 



Hi 



16V2-K 



5vP 



v 



Cp0 2 ,4-i + SffO^A- 



1 - 4/d- + 12k 



(19) 



where V = W or Z, 5y is 2 for the W boson and 1 for the Z boson, j3y is the velocity 

/m 2 H .. Note that 



m 



v 



of the gauge boson in the rest frame of the Higgs boson, and k\ 
two decay widths are of the same form and for m#. ^> mz the ratio of W + W~ to ZZ 
rates approaches 2. Clearly, since the decay widths are proportional to the third power of 
the Higgs boson mass itself, these decay channels are expected to contribute dominantly 
to the decay widths of the heavy Higgs bosons. On the contrary, the partial width of the 
decay Hi — > VV is strongly suppressed or vanishing due to kinematics, thus not playing a 
dominant role even if the decay channel with one off-shell gauge boson is considered. In 
the Born approximation, the intermediate Higgs boson is the pseudoscalar boson so that 
the decay H 2 — > ZZ is not allowed. However, beyond the Born approximation, this mass 
eigenstate will have a large portion of the light or heavy scalar Higgs states and as a result 
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there will be a large enhancement for the branching ratio for non-trivial values of the CP 
phase $ in this specific decay channel unique to the CP-noninvariant theory. 

Thirdly, the heavy Higgs-boson decays into a light Higgs boson and a gauge boson can 
be classified into the charged decays Hi — > H ± W T and the neutral decays Hi — > HjZ. The 
width for the decay if; — > H ± W T is given by 



T{Hi -> ff^T) 



G F m 3 Hi 



cpO 



3,4-i 



4-t 



A 3/2 (l, 



Hi / 



(20) 



where A is the two-body phase space function; X(x, y, z) ■ 
Similarly, the width for the decay Hi — > HjZ is given by 



x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. 



T(Hi -> i^Z) 



G F m 3 Hi 



8^/2 



71 



A 3 / 2 ^1, 



Hi 



mt 



2 ' 2 



(21) 



The branching ratios for the two-body decays H 2 — > H\Z can be sizable in specific regions of 
the SUSY parameter space, especially, for small values of tan (3 and below the tt thresholds. 

Fourthly, the heavy Higgs bosons can decay into a pair of lighter scalars, Hi — ► HjH^ 
and into a charged Higgs-boson pair Hi H + H~ , if they are kinematically allowed. The 
width of the latter decay channel is given by 



H + H- 



16iimH, 



(22) 



with (3jj± the velocity of the charged Higgs boson in the rest frame of the decaying neutral 
Higgs boson Hi and v 2 = v\ + v\. On the other hand, the width for the former decay channel 
Hi -> HjH h is given by 



T(Hi 



HjH k ) 



V \-A-ijk\ ftijk j 



(23) 



where rjjk is 2 for j 7^ k and 1 for j — k, and (5. 



13k 



\ l / 2 (l,m 2 H Jm 2 Ht ,m 2 H Jm 2 Ht ). 
Finally, the decays of the Higgs bosons into two gluons Pp , which occur through quark 
one-loop diagrams may be sizable. As a matter of fact, the two-gluon partial width of the 
Higgs bosons is important because it is the main production mechanism for the Higgs boson 
production at the pp collider LHC. Moreover, since gluons couple to Higgs bosons via heavy 
particle loops, the two-gluon widths are sensitive to heavy particle masses, standard and 
also supersymmetric, well above the Higgs masses themselves. It is a little involved but 
straightforward to calculate the width for the Higgs boson decay into two gluons, which can 
be expressed in the following compact form 



r(fl"i - 99) 



GfOl 2 s w? Hi 
16\/27r 3 



E 



where r = ni 
defined by 



\, /Am 



%\ \F{r) 



R 



RrCl 



|G(r)| s 



(24) 



j, / runs for all colored fermions, and the functions F(t) and G(r) are 
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F(r) 



_-i 



G(r) = r-V(r) , 



(25) 



with 



f(r) = 



arcsin vr 



log 



— ITT 



for t < 1 , 
for r > 1 . 



(26) 



The QCD radiative corrections J2]J to the two-gluon channel may be sizable and are built 
up by the exchange of virtual gluons, gluon radiation from the internal quark loop and the 
splitting of a gluon into two unresolved gluons or quark-ant iquark pair. Although they are 
large enough to nearly double the partial width, we neglect the corrections in our numerical 
analyses because after all the two-gluon channel remains as a sub-dominant decay channel 
for all the neutral Higgs bosons. We note in passing that the two-photon channel and the 
other similar one-loop decay channels have much smaller branching ratios so that they are 
not included in our present discussion. 



C. Numerical results on the branching ratios 

In this section, we present a detailed numerical analysis of the total decay widths of 
each neutral Higgs boson, the branching ratios of each decay mode of the Higgs boson with 
respect to each Higgs-boson mass for various values of the CP phase $. 

Figure 2 shows the total decay widths for the three neutral Higgs bosons with respect to 
each Higgs-boson mass; the left (right) solid line is for the lightest (heaviest) Higgs boson 
and the dashed line for the intermediate Higgs boson. The decay width T(Hi) is almost 
independent of the CP phase, but the other two total widths T(H 23 ) are strongly dependent 
on the CP phase at the lower ends of their allowed masses. The widths are rather small in 
size so that they may not be very hard to be determined at e + e~ colliders, but they may be 
detectable at a fi + pT collider with a very precise energy calibration. 

Figure 3 shows the partial branching ratios of the lightest Higgs boson decays with respect 
to the mass for several values of the CP phase $. Since the lightest Higgs boson mass is 
very strongly restricted to be less than 130 GeV, the main decay channels are into bb, t + t~ , 
gg, and cc. What can be immediately noted in Figure 3 is that the pattern of the branching 
ratios is (almost) independent of the CP phase. We have confirmed quantitatively that this 
feature is not much changed even if the possible off-shell decay channel Hi —>■ W + *W~ is 
included. Numerically, we find that the CP phase $ tends to reduce the contribution from 
the off-shell decay mode. 

On the contrary, as shown in Figure 4 the partial branching ratios for the H 2 decay 
channels with respect to the mass are very sensitive to the CP phase $. The upper 
solid line in each figure frame denotes the sum of the branching ratios of the fermionic decay 
modes. The dot-dashed line is for the decay channel H 2 — > HiZ, two dotted lines for the 
channels H 2 — > W + W~ (upper line) and ZZ (lower line), respectively, and the dashed line 
for the decay channel H 2 — > HiHi. Finally, the lower solid line is for the two-gluon channel. 
Concerning those branching ratios, there are several interesting points worthwhile to be 
mentioned. 
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• Through the whole range of the allowed H 2 mass, the fermionic decay channel remains 
as one of the dominant channels. However, the decay channel H 2 — > H\Z also becomes 
a dominant channel if the channel is open. 

• The decay channel H 2 — > ti overwhelms all other decay channels as soon as it is 
allowed. 

• If the CP phase $ vanishes, the channel H 2 — > H%H\ is prohibited, reflecting that 
the CP-odd state cannot decay into two identical CP-even states when the system 
is CP-invariant. However, the contribution of this decay mode becomes sizable or 
suppressed as the CP phase $ varies from zero. 

• The decay channels H 2 — > W + W~ and H 2 — > ZZ are prohibited for the vanishing CP 
phase. But, they along with the channel H 2 — > H\H\ can be also sizable and become 
the most dominant channel for some values of the CP phase $. 

These interesting features are partly because of a significant mixing between CP-even and 
CP-odd states. As shown clearly in Fig. 5 of the work by Pilaftsis and Wanger |J, the 
couplings of H 2 to R\Z and, in particular, to ZZ increases with the phase $ up to 90° so 
that they clearly enhance the decay channels H 2 —>■ H\Z. On the other hand, the dominance 
of the decay channel H 2 — > H\Hi for a large mass range of the Higgs boson H 2 is mainly 
due to a large enhancement of the coupling strength of the H 2 HiHi vertex with a non- 
trivial phase. This distinct and unique feature may give a strong hint of the existence of a 
non-vanishing CP phase in the stop or sbottom sector. 

When the CP phase is zero, the heaviest Higgs boson is the heavier CP-even state so 
that the decay channels H 3 — > H1H1 and H 3 — > VV (V = W of Z) are among the dominant 
channels. This fact is clearly shown in Figure 5 for the partial branching ratios of various 
H 3 decay channels with respect to the mass m# 3 . In every frame of Figure 5 the upper solid 
line is for the sum of all available fermionic decay channels, two dotted lines for the decay 
channels into two gauge bosons, the dot-dashed line for the decay channel H 3 — > HiZ, the 
dashed line for the decay channel if 3 — > H\H\ and finally the lower solid line for the two- 
gluon mode. These decays of the heaviest Higgs boson H 3 also exhibit several interesting 
features: 

• Before the ti channel is open, the channels H 3 — ^ HxH% and H 3 — > VV are dominant, 
while they are overwhelmed by the ti mode as soon as the latter mode is available. We 
note the similar behavior of B(H 3 — > VV) and B(H 2 — > HiZ) presented in Figure 4, a 
natural consequence from the relation (|12j) among the couplings of the neutral Higgs 
bosons to gauge bosons. 

• The channel H 3 — > H 2 Hi is always prohibited, which reflects that three Higgs mass 
eigenstates are originated from two CP-even states and one CP-odd state, and there 
do not exist any interactions between two different CP-even states and one CP-odd 
state. 

• The channel H 3 — > H 2 Z is not kinematically allowed for the whole mass range even 
though there exists the coupling for the decay mode. 
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• The decay channel H% — > iJiZ is prohibited for the vanishing CP phase, but the 
contribution of the decay channel to the total width becomes sizable as the CP phase 
varies from zero. Moreover, this mode overwhelms the other decay modes near the 
lower tail of the allowed H 3 mass for some nontrivial values of the CP phase. We note 
again the similar behavior of B{H 3 — > H\Z) and BiH.2 — > VV) presented in Figure 4. 

All these interesting features can be understood by noting that the CP-even and CP-odd 
compositions of the heavy Higgs mass eigenstates are strongly dependent on the value of th 
CP phase $. 

To recapitulate, the mass spectrum and the branching fractions of the lightest Higgs 
boson Hi are not so sensitive to the CP-violating phase, but those of the heavy Higgs states 
are very sensitive to the phase. It means that the CP phase may cause in most cases a 
larger mixing of the CP-odd state with the heavy CP-even state than the light CP-even 
state, although the mixing with the light one also may be substantial, in particular for light 
charged Higgs boson masses. From our detailed analysis, it is clear that the CP phase will 
affect the phenomenology of the heavy Higgs-boson decays very significantly so that the 
heavy Higgs boson searches will depend very strongly on the CP phase as well as the other 
SUSY parameters. 



IV. MEASURING CP VIOLATION IN FERMIONIC HIGGS BOSON DECAYS 



In Section |TJ, we have found that the branching ratios of the lightest Higgs boson decays 
are (almost) independent of the CP phase $, and we have noted that for a large value of tan /3 
the neutral Higgs boson masses and the branching ratios of the other neutral Higgs boson 
decays also become insensitive to the CP phase. However, the CP-violation effects caused by 
the loop corrections of soft CP-violating Yukawa couplings can be directly measured by an 
azimuthal CP-odd decay asymmetry in the decays Hi — > t + t~ and/or Hi — > tt [|T^j reflecting 
the correlations of two transversely polarized fermions. In this section, we estimate the size of 
the azimuthal CP-odd asymmetry in the MSSM with explicit CP violation through radiative 
corrections. 

For the sake of discussion, we calculate the helicity amplitude explicitly for the Higgs 
boson decays Hi — > ff by use of the so-called 2-component spinor technique p2| . The 
result is written as 



gm f m Ht 



2mwR 



(T/3fv\ + iRpa 



(27) 



where a and a are for the helicity of the fermions / and /, respectively, with f — rovt. It is 

now clear that the helicity amplitude can be complex when both the scalar and pseudoscalar 

parts co-exist and the interference, signaling CP violation, can be extracted by adjusting the 

polarization of the final fermions. As a result, we obtain the following CP-odd asymmetry 

/. 



in terms of v\, a\, and R 



A 



CP 



21 D__D1 



D 



++\ 



\D. 



2(3 } R f pv)a) 
\P f v)\ 2 + \%a)\ 2 



(28) 
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In order to observe the CP violating spin-correlation of the top quarks or r-leptons, we 
need to look at the correlations among their decay products. The analysis power for the spins 



depends on each decay mode |23|]. For example, if the top quark decays semi-leptonically, 
the spin analysis power is unity, but if it decays hadronically, the analysis power is about 0.4. 
Similarly, if the r decays semi-leptonically, the analysis power is —1/3. And, the inclusive 
hadronic decays of the r is known to give the analysis power of —0.42. Certainly, if one can 
analyze the detailed structure of the hadronic decays, one can obtain better results. 

Although we can do a more comprehensive numerical analysis for evaluating the possi- 
bility of detecting CP violation in the fermion decays, for simplicity we confine ourselves to 
estimating the size of the CP-odd asymmetry in the present work. Since the effectiveness of 
a decay channel in measuring CP violation directly depends on the number of events as well 
as the size of the CP-odd asymmetry, it is useful to define an effective CP-odd asymmetry 
as 



^cp = -4cp V W - ff) ■ ( 29 ) 

With this new effective CP-odd quantity, we can easily estimate the number of Higgs bosons 
required to see the CP-violation effect directly at 1— a level. The number before including 
the detection efficiency and the polarization analysis powers is simply given by 

N H . = — . (30) 
(^cp) 2 



For a more realistic estimate, we refer to the work by Atwood and Soni of Ref. |16] . 



The CP-odd asymmetry is strongly dependent on the value of tan/3. In particular, 
the t + t~ mode is expected to be greatly enhanced as tan/3 increases. In this light, it is 
interesting to consider a large tan/5 case as well, so we consider two values of tan/3 = 3 
and 30 for comparison. The dependence of the CP-odd asymmetry on the charged Higgs 
boson mass m#± for several values of the CP phase in the decays #1,2,3 ~~ * t + t~ is present in 
Figures 6, 7, and 8, respectively, for two values of tan /3; the solid line is for tan (3 = 3 and the 
dashed line for tan f3 = 30. We find that the CP-odd asymmetry in the decay Hi — > t + t~ 
is very sensitive to the charge Higgs boson mass m#± and can be significant near the lower 
tail of its allowed mass range for non-trivial values of the CP phase. Otherwise, the CP- 
odd asymmetry decreases as mjj± increases. This is due to the suppression of the mixing 
matrix element O23 representing the composition of the CP-odd state in the lightest Higgs 
boson. Similarly, near the lower tail of the allowed charged Higgs boson mass, the CP-odd 
asymmetry for H 2 — > r + r~ also is very sensitive to the charged Higgs boson mass. Otherwise, 
the CP-odd asymmetries remain constant instead of decreasing for the whole range of the 
charged Higgs boson mass, but they oscillate between -1 and 1 as the CP phase varies. The 
CP-odd asymmetry for the decay H 3 — > t + t~ is (almost) independent of the charged Higgs 
boson mass, but it is also very sensitive to the CP-phase. On the other hand, the large value 
of the top-quark mass prohibits the lightest Higgs boson from decaying into a top-quark 
pair. So, we consider the decays f/2,3 — > tt, of which the dependence on the charged Higgs 
boson mass m#± for several values of the CP phase is presented in Figures 9 and 10. Near 
the lower tail of the allowed charged Higgs boson mass, the CP-odd asymmetry for H2 — > tt 
is very sensitive to the charged Higgs boson mass as that for H 2 — > t + t~ . Otherwise, the 
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CP-odd asymmetry remains constant for the whole range of the charged Higgs boson mass, 
but it oscillates between -1 and 1 as the CP phase varies, the CP-odd asymmetry for the 
decay H 3 — > ti are almost independent of the charged Higgs boson mass, but it is also 
sensitive to the CP phase. 

V. CONCLUSIONS 

In this paper, we have performed a systematic study of the neutral Higgs boson decays in 
the MSSM with explicit CP violation, which is induced through loop corrections involving 
trilinear CP-violating couplings of the neutral Higgs bosons to scalar top and bottom quarks 
and can be sizable. Certainly, our analysis is far from complete in the sense that any below- 
threshold three-body decay channels and comprehensive QCD corrections, which may be 
significant, are not included. Therefore, a more complete work is needed and it is to be 
reported elsewhere. 

The analysis on the branching ratios of the neutral Higgs boson decays shows that the 
decay pattern of the lightest neutral Higgs boson is almost independent of the CP phase, 
but the decay patterns of two other neutral Higgs bosons are very sensitive to the phase. 
Since the explicit CP violation is realized through the scalar-pseudoscalar mixing between 
the CP-odd Higgs state and, mainly, the heavy CP-even Higgs state, the decay patterns of 
two mass eigenstates are determined by the size of the mixing. Moreover, the dependence 
of each decay mode on the CP phase is different. Therefore, it is naturally expected that 
as the CP phase varies from zero, the branching ratios change a lot. In particular, for the 
intermediate Higgs boson (which is nothing but the CP-odd Higgs state in the CP invariant 
case) the channels H 2 — > H^H^ and H 2 — > VV (V = W or Z), prohibited for the vanishing 
CP phase and so unique to the CP-noninvariant theory, can be sizable and become the most 
dominant channels for some nontrivial values of the CP phase. Furthermore, for the heaviest 
Higgs state (which is the heavier CP-even state in the CP invariant limit) the decay channel 
H 3 — > HiZ prohibited for the vanishing CP phase may overwhelm the other possible decay 
channels near the lower tail of the allowed H 3 mass range for some nontrivial values of the 
CP phase. 

We have also presented a numerical analysis of the CP-odd asymmetry that can be 
extracted through spin correlations of the fermions in the fermionic decays of the neutral 
Higgs bosons, H { — > t + t~ and Hi — > ti. The CP-odd asymmetry is very sensitive to the 
value of tan/5, in particular, for the decay Hi — > t + t~ . For a large value of tan/3, the 
branching ratio for the decay channel is enhanced and moreover the CP-odd asymmetry 
becomes significant for small Higgs boson masses. However, we have found that the CP- 
odd asymmetry for the Higgs boson Hi decreases as the charged Higgs boson mass m H ± 
increases increases due to the suppression of 23 . Except the lower tail of the allowed charged 
Higgs boson mass, the CP-odd asymmetries for both H 23 ~~ * t + t~ and H 2>3 - ► ti remain 
constant for the charged Higgs boson mass, but On the whole, the CP-odd asymmetry can 
be saturated to be 1 so that there is a great chance to see CP violation due to radiative 
corrections directly through the spin correlations of the tau leptons and/ot top-quarks in 
the fermionic decays of any neutral Higgs boson. 

To conclude, the analysis presented in this paper clearly demonstrates that the explicit 



14 



CP violation in the MSSM Higgs sector may modify the Higgs boson decays significantly 
so that it will have a great impact on Higgs searches, and give a great chance to measure 
the effects directly through the CP-odd spin correlations in the decays Hi — > t + t~ and/or 
Hi -> it. 
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APPENDIX 



In this Appendix, we present the explicit form of the coefficients b a and the fully- 
symmetric trilinear Higgs self-couplings (a — 1, 2, 3) in terms of the quartic couplings 
Aj (i — 1 to 8), which are needed for the Higgs boson decays into a charged Higgs-boson 
pair and a neutral Higgs-boson pair, respectively. The coefficients describing the couplings 
of the neutral Higgs boson to a charged Higgs-boson pair in the weak eigenstates are given 
by 

h = 2s ^1 (Age 2 *) - s%X - c 2 X (Are*) , 

b 2 = 2sJ C/? Ai + 4 A 3 - sJ C/3 A 4 - 2s}c f} H (A 5 e 2i «) + S/3 (s 2 - 2c 2 )K (A 6 e*) + s p cpl (\ 7 S) , 

63 = 2c 2 S/3 A 2 + 4A3 - c 2 S/3 A 4 - 2c 2 S/3 ft (A 5 e 2i «) + c^spl (A 6 e*) + C/3 (c 2 - 2s 2 )ft (\ 7 S) , (Al) 

and the symmetric coefficients a Q/ 3 7 describing the self-interactions of three neutral Higgs 
bosons by 

a m = S(S c p l(\ 5 e^) - ±s}l (A 6 e*) - ^1 (\ 7 S) , 
a 112 = i sJc^Ax + ^4(A 3 + A 4 ) - c p {l + sj)ft (A 5 e 2i? ) 
+i^( S 2 - 2c 2 )ft (A 6 e*) + l -cls p K (\ 7 S) 

c 2 S/3 A 2 + -4(A 3 + A 4 ) - s p (l + 4)11 (A 5 e 2 ^) 
+^ 2 c^ (A 6 e^) + i C/3 (c 2 - 2 S 2 )^ (Are«) 
- S/3C/3 J (A 5 e 2 ^) - 1(1 + 2c 2 )J (A 6 e^) , 
-2J (A 5 e 2i «) - C/3S/3 (J (A 6 e l «) + J (A 7 e* 
- S/3 c /3 j(A 5 e 2 ^)-i(l + 2 S 2 )j(A 7 e^) 



aii3 = 

^122 
«123 
a 133 
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= c (i \ 1 + ]-s^l (A 6 e*) 



a 222 = 


c^Ai + 




1 


ri 


^223 — 


3 






1 


ri 


a 233 — 


3 




a 333 — 


s /3^2 + 



3 TZ (A 5 e 2i «) + -s p K (A 7 e 



(A2) 



The remaining 17 coefficients of a aj g 7 can be obtained by the symmetry properties a a/ 3 7 = 
d/3-ia = a-yap = CLa^p = CL-yf3a = CLf3aj- We note that 61, am, Oi22, Qi23, and 0133 come only 
from the imaginary parts of one-loop corrections to the effective Higgs potential so that 
these coefficients vanish in the CP-invariant theory. 
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FIG. 1. Higgs-boson masses with respect to the charged Higgs mass m#± for different values 
of the CP phase $ with the fixed values for other SUSY parameters in Eq. (|l^) . The value of tan (5 
is taken to be 3 and an artificial experimental condition mH 1 > 70 GeV is imposed. The lower 
solid line is for the lightest Higgs boson Hi, the upper solid line for the heaviest Higgs boson H3, 
and the dashed line for the intermediate Higgs boson Hz- For reference, the charged Higgs mass is 
presented with the dot-dashed line. 
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FIG. 2. Total decay widths for the three neutral Higgs bosons with respect to each Higgs-boson 
mass for several values of the CP phase <&. The left (right) solid line is for the lightest (heaviest) 
Higgs boson and the dashed line for the intermediate Higgs boson. 
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FIG. 3. Partial branching ratios of the lightest Higgs boson decays with respect to the mass 
rriH 1 for sevaral values of the CP phase <J>. The upper solid line is for the channel Hi — > bb, the 
dashed line for the channel Hi — > r + r~, the dotted line for the channel Hi — > cc and the lower 
solid line for the two-gluon mode Hi — > 35. 
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FIG. 4. Partial branching ratios for the H2 decays with respect to the mass m# 2 for several 
values of the CP phase <£. The upper solid line in each figure frame is for the sum of the branching 
ratios of the fermionic decay modes. The dot-dashed line is for the decay channel H2 — > H\Z, two 
dotted lines for the channels H2 — ► W + W~ (upper line) and ZZ (lower line), respectively, and the 
dashed line for the decay channel H2 — > H\H\. Finally, the lower solid line is for the two-gluon 
channel. 
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FIG. 5. Partial branching ratios for the H% decays with respect to the mass mu z for several 
values of the CP phase <£. The upper solid line in each figure frame is for the sum of the branching 
ratios of the fermionic decay modes. The dot-dashed line is for the decay channel H3 — > H\Z, two 
dotted lines for the channels H3 — > W + W~ (upper line) and ZZ (lower line), respectively, and the 
dashed line for the decay channel H3 — > H\H\. Finally, the lower solid line is for the two-gluon 
channel. 
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FIG. 6. The CP-odd effective asymmetry A l CP with respect to the charged Higgs boson mass 
m H ± for several values of the CP phase <E> in the decay Hi — > r + r~ . The solid line in each frame 
is for tan/3 = 3 and the dashed line for tan/? = 30. 
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FIG. 7. The CP-odd effective asymmetry A l CP with respect to the charged Higgs boson mass 
m H ± for several values of the CP phase $ in the decay H 2 — > r + r~ . The solid line in each frame 
is for tan/3 = 3 and the dashed line for tan/? = 30. 
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FIG. 8. The CP-odd effective asymmetry A l CP with respect to the charged Higgs boson mass 
m H ± for several values of the CP phase $ in the decay #3 — > r + r~ . The solid line in each frame 
is for tan/3 = 3 and the dashed line for tan/? = 30. 
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FIG. 9. The CP-odd effective asymmetry A l CP with respect to the charged Higgs boson mass 
m H ± for several values of the CP phase $ in the decay H 2 — ► ti. The solid line in each frame is 
for tan (3 = 3 and the dashed line for tan (3 = 30. 
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FIG. 10. The CP-odd effective asymmetry A l CP with respect to the charged Higgs boson mass 
m H ± for several values of the CP phase $ in the decay H 3 — > tt. The solid line in each frame is 
for tan (3 = 3 and the dashed line for tan (3 = 30. 
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